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ABSTRACT 

The investigation of symmetry nonrestoration scenarios has led to a controversy, with cer- 
tain nonperturbative approximation schemes giving indications in sharp disagreement with 
those found within conventional perturbation theory. A Rayleigh-Ritz variational approach 
to the problem, which might be useful in bridging the gap between perturbative and nonper- 
turbative viewpoints, is here proposed. As a first application, this approach is used in the 
investigation of a Z2X .^-invariant thermal field theory with two scalar fields, placing partic- 
ular emphasis on the region of parameter space that has been claimed to support symmetry 
nonrestoration. 
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The subject of temperature-induced phase transitions[l-3] in relativistic quantum field 
theories has been extensively investigated over the last twenty years. In particular, transi- 
tions from a high-temperature symmetric phase to a low-temperature phase in which some 
symmetries are spontaneously broken are a crucial ingredient in most modern cosmological 
scenarios, and have been shown to be realized in large classes of thermal field theories. 

The possibility of symmetry nonrestoration (SNR) at high temperatures[3-7] (or tran- 
sitions from a high-temperature broken-symmetry phase to a low-temperature symmetric 
phase 0) could also have interesting phenomenological implications, most notably allow- 
ing to circumvent the monopole problem in certain Grand Unification Theories; however, 
the investigation of SNR scenarios has led to controversy, with certain nonperturbative 
approximation schemes giving indications in sharp disagreement with those found within 
conventional perturbation theory. Specifically, whereas perturbative analyses find that some 
models, possibly of phenomenological relevance, can support SNR upon appropriate (how- 
ever ad hoc) choice of the available parameters, the corresponding analyses within certain 
nonperturbative approximation schemes @ indicate that symmetry is inevitably restored. 

Recently, some progress has been made toward bridging the gap between perturbative 
and nonperturbative results on SNR. This has been attained via the use of improved pertur- 



bative techniques [|10|, [UJI in which, while preserving the general structure of the perturbative 
expansion, some nonperturbative features of the theory are effectively taken into account. 
The preliminary results of this improved perturbative approaches have indicated [10, 11] that 



the conventional (unimproved) perturbative techniques overestimate the U SNR parameter 
space" (the region of parameter space capable of supporting SNR), and it is reasonable 
to interpret these results as suggesting that further improvements in the accuracy of the 
approximations would ultimately lead to the conclusion that the SNR parameter space is ac- 
tually empty, just as predicted within the nonperturbative approximation schemes adopted 
in Ref.||. This expectation is encouraged by the findings of related studies on the lattice] 12] 



In this Letter, I discuss techniques which can be useful in making further progress in 
the direction proposed in the Refs. JTDj, [TJJ, and, as a first application, which also serves as 
illustrative example, I use them in the investigation of the two-scalar-field theory of Euclidean 
Lagrange density 

L = \(d^)(d^) + \(d^)(d^) + W + W + ^ + ^ _ , (i) 

which is Z2XZ2 invariant [(<& — ► — $) x (^ — ► — \&)], and is among the strongest candidates 
as a model supporting SNR. 

Within conventional perturbation theory the model ([I]) is found to support SNR whenQ 



A(j> < A<j>v]> < y A$A\|/ (2) 

or 



Av]> < A<j>vji < y A(j)Av]/ . (3) 

The lower bound on A$^ in the above @ and @j is set by the requirement that, in the high 
temperature limit, one of the fields have imaginary dressed (thermal) mass, which within 

perturbation theory is found to be given by M ~ J (A$ — A$^)T 2 and Q ~ ^ (A^ — A$*)T 2 

respectively for $ and ^> . (Symmetry is not restored if the square of the dressed mass, i.e. 
the second derivative of the effective potential at the origin, is negative.) The common upper 

1 Note that here and in the following A$ and A* are assumed to be positive, so that there is at least a 
chance of having a stable theory, and all the A's are assumed to be small, so that perturbation theory has 
at least a chance of giving consistent results. 
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bound on A$^ given in (fj) and (0) follows from (naive 
in light of the improved perturbative analyses of Refs. 



TT|) stability analysis 0. However, 
HI it appears that conventional 



perturbation theory underestimates the lower bound of The dNR parameter space, since the 

[□J] that the condition A$ < A$^ (or equivalently 



inclusion of higher order effects shows flO 



is not sufficient to render one of the dressed masses imaginary, and might also 
overestimate the upper bound of the SNR parameter space, since the higher-order analysis 
of ultraviolet structures indicates [|TTJ that the condition A$^ < \/A$A$ is necessary but not 
sufficient for stability. 

While clearly a lot could be learned by further extending/improving the analyses pre- 
sented in Refs. TO, |TTJ, this objective is presented with serious technical hurdles, primarily 
as a result of the necessity to study complicated self-consistent equations for the dressed 
masses (which pick up a momentum dependence beyond leading order). In this Letter, I use 



a variational approximation to estimate the leading correction to the results of Refs. (TO, ITT 



Specifically, I am interested in the way this leading correction modifies the result for the 
lower bound of the SNR parameter space. If a more stringent bound was found, one would 
be presented with a scenario in which each next step of improvement of the analysis appears 
to lead to a smaller SNR parameter space, pointing toward the validity of the above men- 
tioned conjecture that the SNR parameter space is actually empty. However, as reported 
below, I find that instead the leading correction to the result of Refs.JTU], O on the lower 
bound of the SNR parameter space is in favour of SNR. This suggests that additional steps 
of improvement of the analysis might lead to alternating contributions, ultimately resulting 
in a lower bound of the SNR parameter space which does not differ much from the one ob 
tained in Refs. If 
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and therefore is still consistent with SNR in the Z 2 X Z 2 model under 

consideration. 

Since the analysis here presented is purely concerned with the lower bound of the SNR 
parameter space, my result does not affect the possibility that the SNR parameter space 
might be found to be empty once the issues of stability (which are relevant for the upper 
bound of the SNR parameter space) raised in Ref.[|l]J are fully investigated. This however 
appears to require several nontrivial steps, starting with the identification (within a study of 
the type reported in Ref . [|13[1 for the Z 2 model) of the conditions rigorously necessary (and 
sufficient) for the stability of the Z 2 X Z 2 model. This challenging program is left for future 
work. 

In preparation for the derivation of the announced result on the lower bound of the 
SNR parameter space, let me start by revisiting the results of Refs. [ITU], ITTJ. I do this in the 



framework of the CJT formalism of the effective potential for composite operators [11,14- 
16], which is ideally suited for a variational analysis of the type I discuss here. I choose to 
work in the imaginary time formalism of finite temperature field theory, wherein the fourth 
component of momentum is discretized, = iimT (n even for bosons), and I adopt the 
notation 

n=— oo \ "J 

I also exploit the fact that the two SNR scenarios emerging within conventional perturbation 
theory, (fj) and (0), are simply related by renaming of fields and couplings, and therefore, 
for the purpose of the analysis here presented, it is sufficient to consider A's such that 



A$ < A<j>v]> << A\j/ 



(5) 



which allows to investigate the lower bound of the SNR parameter space in the candidate 
SNR scenario @. This is useful since @ implies|Tl[] that all the interesting structures^ of 
the effective potential V(<f>,if}), corresponding to the shifts {$, — > {$ + 0, \I/ + ip}, are 

2 Clearly, the setup (|^) implies that at high temperatures the minimum of the effective potential along 
the ip direction is at %j) = Q. 
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to be found in its projection on the ■0 = axis. Therefore, in the following I concentrate on 
V(d>,if) = Q), i.e. shifts {$, — > {$ + 0, and, for short, I adopt the notation V((f)) for 

In Ref.fllTl this potential V(<fi) was investigated within the "bubble approximation", which 
is the lowest nontrivial order[14-18] of approximation of the effective potential in the CJT 
formalism. The result of Ref.[|n|] can be written, using the (renormalized) gap equations, 
as3 



V 



+ ^ + Q[M] + Q[Q] + + ^0 2 - ^-)P[M] 



24 
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M 2 
Y 



P[M]P[0] , (6) 



where M and f2 are the dressed masses of the fields $ and \l/ respectively, and are determined 
by the gap equations 



M 2 
Q 2 



2 i j 2 

m H <p 

2 ^ 



2 i J.2 

w H — —<p 



A 



P[M] 



(7) 



In the above 



P is the "tadpole" 



X 2 , X 2 
In — 

/i 2 



! + X 2 167T 2 

and Q is the "one loop" @ 

X 4 



(2^ 



k| 2 +X 2 1 — exp 



d 3 k 
(2^ 



In 



1 — exp 



|k| 2 +X 2 




(9) 



(10) 



By stopping at the level of the bubble approximation used in Refs. [jlT] (and effectively 
used in Refs. [pT0]| although the formalism adopted there is different) one ends up ignoring 
contributions to the second derivative with respect to of the effective potential at the 
origin (i.e. the square of the dressed mass of the field $ at = 0) which are of order A 2 In A. 
Taking this into account, and using the following known "high temperature" (small X/T) 
expansions 



rn2 

P[X] ~ — 
L J 12 



XT 

47T 



X 2 , X 2 
In 

16tt 2 T 2 



^ n tt 2 T 4 X 2 T 2 X 3 T X 4 , X 2 

Q\X\ ~ + oln— r , (12) 

^ L J 90 24 12vr 64vr 2 T 2 v ; 



3 Since, for the reasons discussed above, I am considering only the dependence on higher order terms 
of the lower bound of the SNR parameter space, which is quite decoupled Jll| from the issues of stability 
and ultraviolet structure, for brevity my analysis here does not provide details on renormalization (which 
the reader can find in Ref.|ll|), and also ignores the small cut-off dependent contribution to the effective 
potential that is present when the theory, as appropriate for a trivial theory p"5| , is considered with a 
finite cut-off iill. 
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one finds that the Eqs.(^)-(^) lead to the following result for the second derivative of the 
effective potential at the origin 



2 



T 

V"(0) ~ m 2 + (A$-A^) — + A$* v /A*-A<j,*^ 7 =- , (13) 



where indeed terms of order A 2 In A and higher have been dropped, and the relations @ have 

been used in extracting the leading contribution of order AvA. Notice that, in the scenario 
(H), the contribution to ([3~3|) of order A, which can be estimated in conventional (unim- 
proved) perturbation theory, is negative, i.e. works in favour of SNR, but the contribution 

of order \y/\, whose evaluation requires to take into account the bubble diagrams [p 



is 



positive. Most importantly the positivity of the correction of order Xy/\ renders insufficient 
for SNR the condition A$* > A$, which sets the lower bound of the SNR parameter space 



within conventional perturbation theory. This is the basis for the recent observation ||10|, 11 



that the bubble-improved perturbation theory predicts a smaller SNR parameter space than 
conventional perturbation theory. 

The improvement in the determination of the SNR parameter space that I report in this 
Letter relies on an estimate of the leading correction to (|T3"D , which is of order A 2 In A (notice 
that |e| 2 <|e| 2 ln|e|<|e| for small e). I obtain this result via the use of variational techniques 
in the framework of the CJT formalism, which allows to go systematicallyfll, 14-16] beyond 
the bubble approximation. In the CJT formalism the effective potential V is obtained as 
the solution of a variational problem for the effective potential for composite operators W: 



V(4>) = W[cf>;M( ( f>;k)} , (14) 

■ 5W[(f);M(k)} 
. 5M{k) 

A rigorous definition of W can be found in Refs. [14-17], but for the purposes of the present 
Letter it is sufficient to observe that W admits a loop expansion, with Ai(k) appearing as 
the (dressed) mass matrixF]: 



SAIk?i=M(th-k\ 



(15) 



where Vt ree (0) is the tree-level (classical) potential, m tree is the tree- level mass, and W* 
is given by all the two-particle-irreducible jTJJ vacuum-to- vacuum graphs with two or more 
loops in the theory with vertices given by the interaction part of the shifted ($ — > $ + <p) 
Lagrangian and propagator set equal to G(k), with [G _1 (fc)] a b = S a bk 2 + [M. 2 {k)] ab . 

The bubble approximation, which is the lowest nontrivial order[14-18] of approximation of 
the effective potential in the CJT formalism, is obtained by including in W* only the "double- 
bubble diagrams", i.e. diagrams made of two rings touching at one point. Importantly, the 
contributions of double-bubble diagrams have the form 

W y > - whoTTY* = Plx]PlY] • (17) 

4 In the present Letter the formalism is set up in terms of the dressed masses, whereas in Ref . [pj] it was set 
up in terms of the dressed propagators. The two setups are equivalent (if no assumption is made concerning 
the momentum dependence of the dressed objects), but the discussion in terms of the dressed masses guides 
more directly toward the type of Rayleigh-Ritz approximation to be discussed later. 
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which does not involve any flow of momenta from one loop to the other. The possibility 
to make substantial analytic progress, ultimately leading to the result ([]), is a peculiarity 
of the bubble approximation, which in particular involves dressed masses that are exactly 



momentum independent [ 11 , 18 1. As soon as one goes beyond the bubble approximation the 
analysis becomes much more complicated. The dominant higher-order contributions, already 
relevant at the order A 2 In A here under consideration, come from "sunset" diagrams 



Vq{X, Y,Z) jzjz + x2 ^ + y2 ( k+ p y + Z 2 > 



which arise from the contraction of two copies of any three-point vertex present in the shifted 
({$, ^} — > {$ + <ft, Lagrangian. In these diagrams momentum does flow from one loop 
to the other leading to the loss of all simplifications encountered in the bubble approxi- 
mation. In such cases, analytic^ progress requires that the effective potential be evaluated 
approximately. The fact that in the CJT formalism the effective potential is obtained as the 
solution of a variational problem renders available the "machinery" of variational approxima- 
tions. In particular, one can resort to approximations^ of the Rayleigh-Ritz type|T~ 



m 



which, rather than considering arbitrary variations of the dressed masses A4 a b(k) } one takes 
a parameter-dependent ansatz for M. a b{k) and vary only the parameters. This variational 
approach was recently used in R ef.||18|| in a related study of gauge theories at finite temper- 



ature, and has been well received [|1 91, |20||, although it is hard []T8|, [L9| to establish its range of 
validity in the framework of the (thermally troublesome") gauge theories. The non-gauge- 
theoretical model considered in the present Letter is, however, closer to the contexts in which 
Rayleigh-Ritz variational approaches have been traditionally used, and one can expect the 
usual arguments^] for its reliability to hold. Therefore, a Rayleigh-Ritz variational approach 
can be useful in investigations of the issue of restoration/nonrestoration of symmetry in the 
Z 2 x Z 2 model, and, as a first step in that direction, I use it here in estimating the leading 
non-bubble contribution to the effective potential, which comes from the sunset diagrams. 

Before proceeding with this estimate, let me observe that within the bubble approxi- 
mation one could have already resorted to a Rayleigh-Ritz approximation. In fact, as the 
reader can easily check, by substituting in the bubble-approximated CJT effective potential 
for composite operators the following variational ansatz for the dressed masses 

M ab = S al 5 bl M + 6 a2 5 b2 Q , (19) 

one obtains a Rayleigh-Ritz bubble effective potential which, after stationarizing with respect 
to the (dimensionful) parameters M and Q, reproduces Eq.(0) within the accuracy achievable 
in the bubble approximation, with M and f2 identified with the solutions of (H)~(P1)- m f ac t 
this is encoded in the fact that the bubble effective potential obtained in Ref . [|Tl[f is stationary 



(again up to terms which are however not reliably determined by the bubble approximation) 
with respect to variations of M and Q in the neighborhood of the solutions of the gap 
equations (0)-(§). In this sense the Rayleigh-Ritz variational approach is exact within the 

The variational problem in which the CJT formalism casts the evaluation of the effective potential is 
well suited for exact numerical analysis even when non-bubble contributions are included, although, as 
a consequence of triviality, some care is needed in the handling of ultraviolet structures. 

6 The reader should keep in mind that there are two levels of approximation involved in the analysis 
presented in this Letter. To begin with, I am working within truncations/approximations of the effective 
potential that arise within the formalism of the loop expansion of the CJT effective potential. When this loop 
expansion is truncated at the "bubble level" the calculation can be completed without further approximation. 
If one goes beyond the bubble approximation, to the "sunset level" , even the truncated expression of the CJT 
effective potential cannot be analyzed exactly (unless numerical techniques are developed), and a further 
variational approximation is needed to estimate analytically the sunset-truncated CJT effective potential. 

7 Discussion of these topics can be found in Ref. |2l[ . 
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bubble approximation, as it should be expected since the variational ansatz (|T^) is just in 
the form of the exact (bubble-approximated) dressed mass matrix. 

In order to estimate the contribution of the sunset diagrams within the Rayleigh-Ritz 
approach it is sufficient to take the following three steps. First, add the relevant sunset terms 
to the bubble-approximated CJT effective potential. Then, evaluate the resulting expression 
for the effective potential with a physical parameter-dependent ansatz for the dressed masses. 
Finally, stationarize with respect to the parameters characterizing the ansatz. 

Clearly in this higher-order analysis the result of the Rayleigh-Ritz approach is not ex- 
act, but, through the variational procedure, it should provide a reliable plj way to encode 



into a handleable set of parameters the bulk of the effect of the sunset- induced momentum 
dependence of the dressed masses. Actually, as a result of the observation above concerning 
the exactness of the ansatz (|T9| ) at the bubble level, one can expect that the ansatz ([[£]) be 
sufficient also to estimate theleading sunset correction; in fact, as long as the bubble contri- 
butions lead over the sunset ones (i.e. at small values of the coupling constants) the dressed 
masses will be only weakly momentum dependent, so that it should be possible to replace 
them with momentum independent "effective masses" to be determined self-consistently. 

Adopting the ansatz ([19) one finds that the handling of sunset diagrams only requires 
to evaluate Vq(X, Y, Z) with momentum independent X, Y, and Z. This has already been 
done in Refs. [^ , E3f, where, in particular, it was shown that the leading high temperature 
contribution, which is of interest for my analysis, is given by 

T 2 / ST \ 

From Eq. (p!6|) one can easily show that the complete "sunset-approximated" (i.e. obtained 
by adding the contributions from the sunset diagrams to the bubble-approximated effective 
potential) CJT effective potential can be written as 



V + ^ + Q[M] + Q[n] + (^ + ^ - ^-)P[M] 



+(y + *f <P 2 ~ %)PP\ + y (P[M]f + ^(PMr - ±f P[M)P[n] 

X 2 X 2 
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v q (m, m, m) - -^(f) 2 v Q (n, n, m) . (21) 



As discussed above, the ordinary effective potential is obtained from this formula by sta- 
tionarizing with respect to variations of M and Q. Then, also using the high temperature 
expansions (pTTj) , (|T2"D, and (pC|) , with some simple algebra one can show that the sought 



contributionjjof order A 2 In A to the second derivative of the effective potential at the origin 
is given by 

Arf 2 , A$M 2 M t \UT\ m + M. A|T 2 

ln (^) - -TT^JT H^f) + -TTT^- M — — — ) + -—r ln(-) ~ 



16tt 2 y T J 16vr 2 y T J 16n 2 v 3T ' 48tt 2 v T 

A$*fi 2 X^X^T 2 24 

— ^-m — ) ~ —In — < , 22) 

16vr 2 K T J 768tt 2 A* V ; 

where the indicated approximations hold in the scenario (|j) presently under consideration, 
which also implies M « Q. 

8 Tcrms of order A 2 In A are already present in the bubble approximation (these terms are among those 
listed in Eq.(p2[)), but could not be included JTo[ i n the bubble approximation of the second derivative of the 
effective potential at the origin, given in Eq.(|l3|), since they are not reliably determined when the sunset 
contributions are not taken into account. 
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Eq.(|22|) shows that, as anticipated in the opening of this Letter, the leading 0(A 2 lnA) 
correction to Eq. (pT3|) is negative, i.e. works in favour of SNR. This is certainly of encourage- 
ment for the possibility of SNR, with the proviso about stability mentioned in the opening 
of this Letter. It should also be noticed that the SNR-favouring correction here identified is 
quite small, especially as a result of the numerical prefactor, and, while possibly sufficient to 
help the case of SNR in the toy model here considered, it might be negligible in the context 
of models rel evan t for particle physics, where large gauge-coupling contributions render SNR 
very unlikely |TD]. The Rayleigh-Ritz variational approach here advocated should be useful 
in further investigation of these issues. 



It is a pleasure to acknowledge conversations with G. Bimonte, R. Jackiw, G. Lozano, O. 
Philipsen, S.-Y. Pi, and S. Sarkar. 
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